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A spin and phonon excitations spectrum are studied using quantum Monte Carlo method in 
antiferromagnetic chain with spins S = 1/2 coupled nonadiabaticity with acoustic phonons . 
It is found the critical coupling exists to open gap in the triplet excitation spectrum for any 
phonon velocity. The phase boundaries of delocalized phonons and propagated the bound 
states of magnon and a phonon are calculated. It is shown that the spherical symmetry of 
the spin-spin correlation functions is broken . The magnetic and optical properties CuGeO^ 
are explained without using spin-Peierls transition. 



PACS: 75.10.Jm, 75.30.Ds, 67.20. +k, 63.20.Ry 

The existence and stability of multiquanta bound states are now established for a wide variety 
models with prescribed nonlinearity and exhibited some distinctive observable signatures in 
terms of spatiotemporal correlations. The typical origin of effective nonlinearity in quantum 
system is the coupling of two or more fields. Adiabatic slaving of fields usually results in 
nonlinear Schrodinger models. Realistically nonadiabatic effects must be considered and the 
influences of nonlinearity and nonadiabaticity are inevitably interrelated. Here we consider 
spin-phonon coupled model frequently used to describe the opening of a singlet-triplet spin 
gap and an isotropic drop of the magnetic susceptibility below critical temperature. As a rule 
a optical phonon interacts with spin system and lead to lattice instability. The interaction 
between acoustic phonons and spins has been studied in two limit cases. A weak coupling 
lead to a renomalization of the phonon mode and give rise to four spin interaction. As a 
result of a gap is opened but the other well-defined magnetic and phonon bound states has 
been loosed. A strong coupling dos't correspond to real compounds and the results obtained 
in the intermediate coupling are incorrect. 

A magnetic, elastic and optical properties in the inorganic quasi-one-dimensional CuGeO^ 
compound fail to describe in terms of one model. The thorough structure measurements [jll 
and Cu-nuclear quadrupole resonance experiments have not been found a lattice dimer- 
ization and a softening of optical phonons. The static distortions below critical temperature 
Tc = lAK cannot be ascribed to any pseudomode in the high-symmetry phase 0. The mag- 
neto induced structural effects depend on the used frequency , so strong lattice fluctuations 
have been observed by electron diffractions up to T ~ lOOi^ by optical methods up to 
T ~ 200i^ 0. A frustration interaction due to a asymmetric or a next nearest- neighbors 
antiferromagnetic exchange can cause a gap in the triplet excitation spectrum though no 
experimental proves have existence of these exchanges. The physical nature of the critical 
temperature existence and opening gap in triplet excitation spectrum, the Raman spectrum, 
thermal conductivity in magnetic field remains unclear. 

Numerical methods, such as exact diagonalization and density-matrix renormalization group, 
face the potential difficulty in dealing with a very large Hilbert space of the phonons. Here , 
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we use Monte Carlo approaches restricted to finite chains L = 100, 200 but without any adi- 
abatic approximation and the truncation of the infinite phonon Hilbert space. The method 
[^] is based on a path-integral representation for discrete system in which we work directly 
in the Euclidean time continuum. All the configuration update procedures contain no small 
parameters . Being based on local updates only, it allows to work with the grand canonical 
ensemble and to calculate any dynamical correlation function, expectation values. 
We consider a model Hamiltonian of an spin-phonon system: 

H = j2[J + a{u,-u,+^)][S:SU^ + iStSr^, + SrSt^,)/2]+Mu^/2 + K{u,-u,+,y/2. (1) 

Here S'^'^ are a spin operator components associated with the site i, J > is the usual 
antiferromagnetic exchange integral, a is the spin-phonon coupling constant, Ui is the dis- 
placement in the z- direction, M is the mass of the ion and K the spring constant. Using 
the quantum representation for phonon operators b,b'^, the Hamiltonian becomes 

H = EE[^+«\/^ysi^cos(g(^ - 0mK+bt)][S!SUMStSr^, + SrSt,,)/2]+ (2) 

Y,hn{q)bpg, n{q) = 2^sm{^), uJo = 2^, g = 27m/L, n = 1, 2, ...L 

The phonon frequency , spin-phonon coupling, energy, temperature are normalized on the 
exchange J. The temperature used in calculation is /3 = J/T = 25. Our system consists 
of the two subsystem interacted. The elastic subsystem is described by phonons with the 
number of occupation Uph = 0, 1, 2.. and magnetic subsystem is characterized by the number 
occupation rim = 0,1 and Pauli operators a, a"*" which coincide with S"^ spin operators. We 
start with the standard Green function of the phonon in the momentum q- imarginary-time 
r representation: 

G{q,r) = E I < ^l&^l^ac > pexp [-{E,{q) - Eo)t], (3) 

where |z/ > is a complete set of the Hamiltonian H in the sector of given q , if|z/(g) >= 
Ey{(i)\v{(l) >, H\vac >= Eo\vac >,Eq = 0. Rewriting Eq.(3) as 

G{q,T)= dcuA(g,tu)exp-(cur),A(g,cu) = < z/|6+|mc> |2 (4) 

Jo I, 

one defines the spectral function A{q,uj) . We calculate a one-particle spin Green function 

/•oo 

< a(r)a^(0) >= / duopt{Lo) exY> —{uot), (5) 



where pt{^^) -spectral density function , and two-particle spin Green function associated with 
singlet excitation AS"^ = 

poo 

< ao(r)a]^(r)ai(0)aQ (0) >= / (ico'ps(ci;) exp — (cjr). (6) 

Jo 

Correlation functions are determined on the basis of a complete set of eigenstates of the 
Hamiltonian: 

<0>= ^^<^^\0\^^> ^ (7) 
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where O = {aoaf + aoCti), {bq + 6;j")S'oS'f, {aoaf + aoCti)(&y + &^)- The longitudinal spin- 
spin correlation function < SfSf^^^ >, the phonon density- density < nph{q)nph{q + p) >, 
the distribution of phonons number nph{q) and magnons number njn{k) as a function of 
momentum are simulated. 

Energies excitation of the two-particles S'^S^_p,p — tt and one- particle excitation differ 
by less than ~ 1% in the isotropic antiferromagnetic chains. The interaction with the 
elastic system lead to qualitatively various dispersion relation. The two-particle excitation 
spectrum became gapped and one-particle excitation gapless for a < ad- The gap energy 
determined from spectral density well fit on the linear dependence As{a) ^ A < J . 
For small spin-phonon coupling a < ccd the phonons are created near the upper bound of 
triplet excitation band Wt with momentum q = 2arcsin^. It follow from the distribution 
of the phonon number versus momentum and the spectral density. The spins number in the 
structure of the spin-phonon quasiparticle h'^_f.^_j^^_ ___^^S].^S].^...S].^,'^ = increase and 
the correlation function between spin and phonon 




= \It7< I COS {q{i - 0.5))(6+ + bq)S!SU, >, (8) 

qualitatively confirms it as illustrated in Fig. la . For a < «c2 the correlation function is 
parameterized as Rbs — ,ujq > ?> Rbs — ,ujo < 3. Correlation function between 
the longitudinal acoustic phonon and the transverse spin components < (6+ -|- 6)ao oi > 
is less by a factor 3-5 than Rbs- Number of two-particle singlet excitation falls as to an 
exponent Ns ~ |exp— f^^]. The formation of the bound states and nonlinear excitations 
causes to break the spherical symmetry of the spin-spin correlation function. The ratio of 
the correlation functions of transverse spin components to longitudinal is parameterized as 
the linear dependence on the spin-phonon coupling 

r< SftST' > ^, 5,2a 

|_^_21.^,„<.„„ (9) 

as presented in Fig. lb. The magnetic susceptibility is represented by spin-spin correlation 
function when the long-range order is absent Xx(y)IXzz = J2r SoS~/ X^r SqS^. Conventionally 
the physical nature of anisotropy result from spin-orbital interaction for spin S = 1/2. 
We suggest a new mechanism the dynamical interaction between the longitudinal acoustic 
phonon and spin. It is important for S- ions. 

For a > Oci the gap in the triplet excitation spectrum simulated from the spectral density 
function is open . The functions p arc plotted in Fig. 2. The gaps energies fit satisfactory 
on the straight line A ~ 0.8(9) (a — ad) in the range of ad < a < ac2- Here exists 
the bound states of spin-phonon quasiparticls similar to bipolaron with the symbolic kind 

Y[? Sf >. The correlation function of the phonon number < 

Nph{q)Nph{q + p) >, presented in Fig. 3b is nonzero at the certain wave number p for a ~ ad 
and for all p at the condition a > ac2- The magnon distribution function demonstrates a 
small oscillations approximately within a five pcrccnts as compared to n„((g = 0,7r) ~ 0.22 
for a > ad and the number phonon distribution function reveals four local maxima for 
a > ac2 and one sharp maximum for a > ctcs. Analysis of the spin- spin correlation function 
plotted in Fig. 3 lead to two important length scales: 

I I ^ / / r COS (Qr) COS (qr) 

I < S^S^ > I ~ -exp(-r/0, r < rc, < S^S^ >~ (-1)" r > Vc- (10) 
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where ^ correlation radius, Q,q,6 fitted parameters. R n/Q may be interpreted as the 
typical size of a cloud of spins and phonons ; ^ may be considered as the width of the 
transition region separating the nearest spin-phonon quasiparticles. These coherent nonlinear 
excitations cause the local extremums of the spectral density phonons and spin excitations 
at the same energies including the spin gap range as shown in Fig. 2. It is possibly a standing 
vibration with the spin kink in sites exists here. The disjoint vibrations are performed at 
the condition that the wave length is changed in two times //a = 2; 4; 8; 16; ... c k = n/l, E = 
wosin{k/2). These estimates of energy are in good agreement with Monte Carlo results up 
to n = 4 because the temperature fluctuations (T = 0.04),,^ < l/a = 64 cut the kinks 
interaction radius. The average number of phonons increase vs. a according to the power 
law 

Na. = 0,0012(— < ^ o,0006(— > t;^, (11) 

here Vph and Vm is the velocity of phonon and magnon. The spectral density of phonons and 
the distribution of phonon number became continuum at a = ac2- The value ac2 is similar to 
critical concentration when dressed phonons effectively "percolate " along the three low-lying 
bands. The correlation function < Nph{q)Nph{q+p) > is not equal to zero for all momentums 
as shown in Fig.3. The transition from one state with localized phonons to delocalized is 
like to Anderson's transition in disordered systems. The width of band increases versus a 
as Wph{ac2) /Wph{aci) ~ 1.5. When the spin-phonon coupling is compared to the energy 
of the boundary of top band of phonon and triplet spin excitation a > ac3 ~ O.Swq the 
soliton lattice is formed with the wave vector of structure being in the range of vr < Q < 27r 
. A nonlianear excitation ,for example, the spin breather is gapless and the spectral density 
of phonons and spin excitations are comparable values at the low energies . The gap in 
one-particle spectrum excitation pass into quasigap. The area of parameters in the plane 
the energy of top band of acoustic phonon uq - the value of the spin-phonon coupling a , 
where a gap in the one-particle excitation spectrum is arisen, is limited by two fitting lines 
ctci — 4^ and ac3 ~ O.Swq ■ The complete filling of the three acoustic bands occurs at the 



~ ^0 



«c2 — 2Wj- 

Now we will apply our results for interpretation and prediction of new effects for the one- 
dimensional spin system CuGeOs with intrachain exchange Jce = 120i^, Ji/Jce = 0.1 0. 
The crystal structure CuSiOs is isostructural to CuGeO^. However CuSiO^ is antiferromag- 
net with Jsi = 21K, Ji/Jsi = 0.14 (Ji— interchain exchange )@|. The interaction between 
spins and acoustic phonon mode in chain leads to gap in the triplet excitation spectrum 
CuGeO^. For the same value of the spin - phonon coupling is [ojq/ J]si » [^jJo/J]Ge and 
(^ci,si > ctci^Ge- The magnon velocity CuGeOs is Vm ~ 1300m/ s, the phonon velocity depends 
on the ultrasound frequency and we take the average value Vph ~ 5000m/s = ^ = 
3,84. The anisotropic elastic constants are cn = 64GPa,C22 = 37.6G'Pa, C33 = 317.3GPa 
[ p!0|] and confirm validity using one dimensional model to study the dynamic of the interacting 
magnetic and elastic systems. According to our results the spin-phonon coupling in CuGeO^ 
has been established a ~ 2, 2 for A/ J = 0.2, a;o = 12. The estimated energies of the bound 
spin-phonon excitation are Egph = 95cm~^, 48cm~^, 24.4cm~^, 12.5cm~^, 6.2cm~^, 3cm~^ and 
are in general agreement with Raman spectrum = 98cm~^, 48cm~^, 30cm~^, 13.2cm~^ 
0. The localized phonon excitation in the second E2^'~^ = 0.22eV and in the third 
E^^^' = 0.34eV^ bands agree satisfactory with data of the optical absorbtion £^2^ = 0.21eV 
[|1^, E^^ = 0,36eV and make clear the physical origin of these peculiarities . Up 



to present time the existence of isolated energy level E^^ ~ 6cm below the triplet gap 
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has been remained unclear. This gap is found by the inelastic neutron scattering 



and the submillimiter resonance |T^. We predict a new resonances at the frequencies 



uo ~ IdiQGHz] ~ 65GHz at low temperatures T < IK. 

The bound spin-phonon excitations will result in additional phonon scattering which depen- 
dents on magnetic field. The energy of spin-phonon quasiparticle normalized on gap value 
Ef^^ / IS. = 0.37; 1.5 is well agreement with temperatures normalized on critical tempera- 
ture Tc = lAK ,T^^ /Tc ~ 0.39; 1.57 related to local maxima of thermal conductivity p6| . 



This effect result from interaction between elastic and magnetic system and magnetostric- 
tion parameter exhibits also two maximum at the same temperatures []T^. The relationship 

of transverse spin-spin correlation function to longitudinal [^^fi^ ^2] ~ 0,9 is associ- 
ated with the anisotropy exchange (Jb/Jc)^^^' = 1,45 which agrees with experimental data 



Jb/Jc = ^b/^c = 1,38 |T^. The sharp rise of g— factor near the transition T > Tg may 
be explained in term of formation of the coherent state of spin-phonon quasiparticles which 
may be broken by thermal phonons. The average number of phonons resulted in spin-phonon 
interaction < Nq >= 1,7 ■ 10~^ in CuGeO^ is equal to number of thermal phonons define in 



terms of Debye approximation with Debye temperature G = 330i^ |jT9| at the temperature 
T = 18ir. 

In conclusion, we summarize the main results. The interaction between magnetic and elastic 
system lead to gap in one-particle and two-particle excitation spectrum. The gapped triplet 
spectrum arises at the certain critical value of the spin-phonon coupling. Bound spin-phonon 
excitation localize for a < and causes a pecuharities set of spectral density of phonon 
and spin excitations. Spherical symmetry between transverse and longitudinal spin-spin 
correlation function is broken. The triplet gap in CuGeO^ arises from the interaction between 
spin and acoustic phonon and this compound has not possessed spin-Peierls ordering. From 
our results the magnetic , elastic and optic properties CuGeO^ is explained. 

The author are grateful to V.A. Kashurnikov, A.S. Mishchenko, N.B. Prokof'ev, B.S. Svis- 
tunov for useful assistance and discussions of Monte Carlo method in the continuous time. 
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Captions to Aplesnin paper S=l/2 " Quantum spin liquid in antiferromagnetic chain S=l/2 
with Acoustic Phonons" 

Fig.l The correlation function of the nearest- neighbor spin-spin and displacement < Aqi{b'^ + 
bq)SfSf^i >, where Agi — y^y'sin | cos {q{i — 0.5)) (a) and ratio transverse spin-spin cor- 
relation function to longitudinal (6) with loq — 10(1), 8(2), 6(3), 1(4) versus spin-phonon 
coupling . 

Fig. 2 The spectral density onc-particle excitations with uq = 8,a = 0.4(1), 2.6(2), 3. 2(3)(a), 
a = 1.6(1), 2(2)(6); and two-particle singlet excitations with ujq = 8,a = 2.6(2), 3.2(l)(c) 
and phonon excitations with loq — S,a — 2(1), 3.2{2){d) . 

Fig. 3 Longitudinal spin-spin correlation function vr. r witha;o = 8, a = 0.8(1), 2(2), 3. 2(3) (a). 
Correlation function of the phonon number normalized on its maximum value versus mo- 
mentum with cuo = 8, a = 1.4(1), 2(2), 2.6(3) (6). 
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